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CLAIMS 

1. Method of determining the velocity V and 
anellipticity r| parameters for processing seismic traces 
in a common midpoint (CMP) gather including an anelliptic 
NMO correction, comprising: 

- a preliminary step to define a plurality of nodes 
(dtn. To) , the said nodes being indicative of parameters 
dtn and Tq representing the NMO correction' for the 
maximum offset and the zero offset travel time in 
hyperbolic coordinates, the said preliminary step being 
followed by 

- for each node (dtn, Tq) defined in the preliminary 
step, the following steps: 

- for static NMO correction of traces in the CMP 
gather as a function of the values of the said parameters 
dtn and to at the node considered., and 

- for calculating the.. semblance function 
associated with the said NMO correction for the. node ' 
considered; and V ' . - 

- for each picked'.. . "time to, a ' step including 
determination of ' the" maximum ' semblance node (dtn (to)', 
T:o(to)) , 

- and a final step to convert the dtn (to) and To(to) 
parameters so as to. obtain the velocity V(to). and 
anellepticity IT (to) laws. 

2. Method according to claim 1, wherein- the nodes 
are defined during the preliminary step in an analysis 
volume (dtn, Tq, to) determined by minimum and maximum 
values respectively [dtn,,i„, dtn^^ax] [Tomin. Xotnaxl and [tomin. 
tomax] of the dtn, Tq and to parametei's - 
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3. Method according to claim 2, wherein, during the 
preliminary step, a corridor [dtn™i„( to) , dtn_(to)], 
tTomin(to), w(to)] for Changing dtn and tp parameters is 
delimited inside the analysis volume as a function of 
plausible velocity V and anellipticity ^ values, the 
nodes (dtn, to) defined for applying the NMO correction 
being then located along the corridor thus delimited. 

4. Method according to any one of the preceding 
claims, further comprising, for each node (dtn. to), a 
stacking step of the corrected seismic traces, following 
the semblance function calculation step. 

5. Method according to claim 4, wherein the stacking 
of corrected traces is done using only near offset 
traces. 

6. Method according to any one of claims 4 or 5, 
further comprising for each picked time, and following 
the step for determining the maximum semblance node, a 
step of checking that values . dtn and Xo ' of the maximum 
semblance node correspond' to a stacking extreme valiie. f or 
the same values dtn and Tq . 

7. Method according' to any one of the' preceding 
claims, further ■comprising, a step of ^ •selectin^^.', and 
adjusting the pickings obtained. following the step 
implemented for determining the maximum semblance node 
(dtn(to). Xo(tc,)) for each picked time to, before the 
conversion step. 

8. Method according to claim 7, wherein the said 
step of selecting and adjusting the pickings comprises a 
step of only retaining pickings dtn and to for which time 
to the highest semblance pickings is greater, than a 
predefined value. 
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9. Method according to claim 8, wherein the said 
step of selecting and adjusting the pickings also 
comprises a step for adjusting the retained pickings dtn 
and To by parabolic interpolations using values about the 
said picked values. 

10. Method according to claim 9, wherein the said 
step of selecting and adjust ingt pickings also comprises 
a step of eliminating retained and adjusted pickings dtn 
and To when it is impossible to calculate the Dix interval 
velocities between the picking considered and higher 
semblance pickings. 

11. Method according to any one of claims 1 to 10, 
wherein the processing applied to seismic traces is an 
NMO correction process implementing a static correction 

15 CORRnmo. 

12. Method according to claim 11, wherein, during 
the preliminary step, the NMO corrections CORR^*, are 
calculated for all nodes (dtn, to) included in the 
analysis volume and all offsets of processed seismic 

2 0 traces. 

13. Method according to claim 12, wherein the NMO 
correction carried out for each node' (dtn, to) consists of 
applying NMO corrections CORR^ calculated during the 
preliminary step. 

14. Method according to any one of claims 11 to 13, 
wherein for a given (dtn, Tq) pair, the static NMO 
correction CORRn„o of a seismic trace with offset x is 
carried out according to the following equation: 

CORR«„o(x)= - TO . L^ dtn(dtnH-2T,) ^, ^^.^^ 
30 represents the maximum Offset in the CMP gather. 
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15. Method according to any one of claims 1 to 10, 
wherein the processing applied to seismic traces is a 
PSTM migration using a static NMO correction CORRps™. 

16. Method according to claim 15, wherein, during 
the preliminary step, ' the NMO corrections CORRpstm are 
calculated for all nodes (dtn and xo) included in the 
analysis volume and all migration offsets inside the 
migration aperture. 

17. Method according to claim 16, wherein the NMO 
correction step carried out for each node (dtn and tq) 
comprises, for each offset class, application of the said 
NMO corrections CORRpstm, calculated during the 
preliminary step on all midpoints inside the migration 
aperture . 

18. Method according to claim 17, wherein the NMO 
correction step carried out for each node (dtn and xo) 
comprises, for each offset class, the stack of the 
corrected midpoints following application of the said NMO 
corrections CORRpstm- ' ' ■• 

19. Method according to any one of. claims 15 to 18, 
wherein, for a given pair (dtn and Xo) , the ' static NMO 
correction CORRpstm is carried out according to the 
following equation: 

CorrpsTM(x) = -Tq + /r,- _ dtn(di,t + 2r,Xx- x + h)' + \^ ^ dtn{din + It A(x - X - hY 

h — 

25 where: .... 

- Xn, represents the. coordinates of the midpoints, 

- X - x„ represents the migration aperture PSTM, 

- h is the half source - receiver offset, 

- xmax is the maximum offset and aperture of the 
30 migi-ation. 
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20. Method according to any one of claims 14 or 19, 
wherein, during the final conversion step, the parameters 
dtn (to) and (tO) are converted" to the velocity law V{tO) 
according to the following equation; 



V = 



dtn{dtt7 + 2To)^ 



21. Method according to one of claims 14 or 19, 
wherein, during the final conversion step, the parameter 
TO (to) is converted to the anellepticity ri(tO) law 

according to n = 

22. Method according to claims 20 and 21, wherein 
parameter dtn is defined with respect to the velocity V 
and anellepticity rj according to the following equation: 




dtn =^^* . /I ^0 



^2 

X 



l + SiiJ (1 + 8ti)V^ 



23. Method according to claim 21, wherein parameter 
TO is defined according . to anellepticity r| accoi-ding to 
the following equation: 

1 + 8T1 

24. Method of characterising a velocity field for 
processing seismic data using a gather of seismic t-races • 
at common midpoint, wherein, for each travel time to for 
a zero offset, a set of parameters dtn and tO is defined, 
representing the NMO correction for maximum offset, and 
the zero offset travel time respectively, in hyperbolic 
coordinates. 



